clear
use pw_coalitions.dta

fp <pceyr>: probit failure maxshift losercinctvc shiftloser lognum loglength <pceyr>, cluster(warnum)

hist maxshift if e(sample), freq
graph save Graph maxshift.gph, replace

hist losercinctvc if e(sample), freq
graph save Graph cincloser.gph, replace

graph combine maxshift.gph cincloser.gph, saving(histograms.gph, replace)

clear
